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ABSTRACT 

Frequencies of intermediate-degree f-modes of the Sun seem to indicate that the solar 
radius is smaller than what is normally used in constructing solar models. We inves- 
tigate the possible consequences of an error in radius on results for solar structure 
obtained using helioseismic inversions. It is shown that solar sound speed will be over- 
estimated if oscillation frequencies are inverted using reference models with a larger 
radius. Using solar models with radius of 695.78 Mm and new data sets, the base of 
the solar convection zone is estimated to be at radial distance of 0.7135 ± 0.0005 of 
the solar radius. The helium abundance in the convection zone as determined using 
models with OPAL equation of state is 0.248 ±0.001, where the errors reflect the esti- 
mated systematic errors in the calculation, the statistical errors being much smaller. 
Assuming that the OPAL opacities used in the construction of the solar models are 
correct, the surface Z/X is estimated to be 0.0245 ± 0.0006. 
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1 INTRODUCTION 

Helioseismology has proved to be a remarkable tool that 
can be used to study the solar interior. Helioseismic inver- 
sion techniques have been used to probe details of the solar 
structure (cf., Gough et al. 1996; Basu et al. 1996; Koso- 
vichev et al. 1997), abundances (cf., Basu & Antia 1995; 
Kosovichev 1996; Antia & Chitre 1997a), and rotation rate 
(cf., Thompson et al. 1996). Most inversions done so far as- 
sume that the radius of the Sun is known precisely and the 
value used is 695.99 Mm (Allen 1976). 

Recent data from the Global Oscillation Network Group 
(GONG), and from the SOI/MDI instrument on board 
the SOHO spacecraft, allow the determination of the 
intermediate-degree solar f-mode frequencies very precisely. 
It has been shown by Antia (1997) and Schou et al. (1997) 
that the ratio of the newly determined solar f-mode frequen- 
cies and frequencies of the f-modes of standard solar mod- 
els is very nearly a constant which differs significantly from 
unity. Since the f-mode frequencies are largely insensitive to 
the detailed structure of the Sun and depend mainly on the 
global parameters, the simplest explanation of this difference 
is that the radius assumed for the models is not the radius 
of the Sun. Antia (1997) estimated the solar radius to be 
695.78 Mm using the GONG data, while Schou et al. (1997) 
using MDI data show that the apparent photospheric solar 
radius (695.99 Mm) used to calibrate the models should be 
reduced by approximately 0.3 Mm. This difference between 
the radius of the Sun and of solar models used as reference 
models for helioseismic inversions can lead to errors in our 



estimation of the solar structure. 

The exact value of the solar radius as determined by 
the f-modes depends to some extent on the surface physics 
used in the models (cf., Antia & Chitre 1997b). Although 
the frequencies of solar f-modes are remarkably insensitive 
to the structure of the solar interior, they do depend to a 
small extent on the details of the surface layers. The modes 
are affected particularly by the formulation used to calculate 
the convective flux and atmospheric opacities. Models con- 
structed with the Canuto-Mazzitelli (1991) formulation for 
calculating the convective flux show that a smaller reduc- 
tion in radius is needed for the model and observed f-mode 
frequencies to match, compared with the reduction needed 
for models constructed with the mixing length formalism. 
This difference is most probably responsible for the appar- 
ent difference of results found by Schou et al. (1997) using 
the data from the SOI/MDI instrument on board the SOHO 
spacecraft, and by Antia (1997) using data from GONG. As 
a result, there is still some uncertainty in the estimate of 
solar radius. However, if models with the same physics are 
used, then the GONG and MDI data give very similar results 
(Antia & Chitre 1997b). 

An inversion for solar structure (e.g., Dziembowski, 
Pamyatnykh & Sienkiewicz 1990; Dappen et al. 1991; Antia 
& Basu 1994a; Dziembowski et al. 1994) generally proceeds 
through a linearisation of the equations of stellar oscillations 
around a known reference model. The differences between 
the structure of the Sun and the reference model are then 
related to the differences in the frequencies of the Sun and 
the model by kernels. Non-adiabatic effects and other er- 
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rors in modelling the surface layers give rise to frequency 
shifts (Cox & Kidman 1984; Balmforth 1992) which are not 
accounted for by the variational principle. In the absence 
of any reliable formulation, these effects have been taken 
into account in an ad hoc manner by including an arbitrary 
function of frequency in the variational formulation (Dziem- 
bowski et al. 1990). Thus the fractional change in frequency 
of a mode can be expressed in terms of fractional changes 
in the structure of the model and a surface term. 

When the oscillation equation is linearised — under 
the assumption of hydrostatic equilibrium — the fractional 
change in the frequency can be related to the fractional 
changes in two of the functions that define the structure 
of the models. Thus, 
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(cf., Dziembowski et al. 1990). Here 8u)i is the difference in 
the frequency u>i of the ith mode between the solar data 
and a reference model. The functions /i and fa are an ap- 
propriate pair of functions like sound speed and density, or 
density and adiabatic index Fi etc. The kernels K[ 2 and 
K\ x are known functions of the reference model which re- 
late the changes in frequency to the changes in f\ and f 2 
respectively; and Ei is the inertia of the mode, normalized by 
the photospheric amplitude of the displacement. The term 
-fsurf results from the near-surface errors. 

The kernels K\ 2 and K 2: i are normally calculated as- 
suming that there is no difference between the radius of the 
Sun and the reference model (see Antia & Basu (1994a) 
for the steps involved in calculating the kernels for the pair 
(c 2 ,p) ). Thus inversions of Eq. (1) are likely to introduce 
systematic errors if the radius of the reference model is dif- 
ferent from that of the Sun or test model. This can be seen 
in Fig. 1 for inversion between two models. Hence, our re- 
sults on solar sound speed are also likely to be modified if 
the radius of the solar model is different from that of the 
Sun, and initial results suggest that this is indeed the case 
(Antia 1997). 

In this work we show how a mismatch of radius af- 
fects solar sound-speed results. Further, since most of our 
knowledge of the solar interior, e.g., depth of the convec- 
tion zone (CZ), overshoot below the solar convection zone, 
helium abundance in the solar envelope, etc., are derived 
from the sound speed, we re-calibrate these quantities in 
view of possible changes in the value of solar radius. While 
determining the helium abundance we have looked into dif- 
ferent sources of systematic errors, like opacities and depth 
of CZ, which had not been adequately examined in earlier 
studies. We also attempt to put limits on the heavy-element 
abundance in the solar envelope. Unlike most previous works 
which utilized oscillation data from the Big Bear Solar Ob- 
servatory, in this work more precise data from GONG and 
SOI/MDI have also been used to determine the depth of the 
solar convection zone and the helium abundance. 

The rest of the paper is organized as follows. In Section 
2, we discuss the change in deduced sound-speed and den- 
sity profiles in the Sun that results from using models that 
correspond to the reduced radius as determined from the 
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Fig. 1 The result of sound-speed inversions between different 
models. The continuous line is the exact sound-speed difference 
between the two models, the points are the results obtained by 
inverting the frequency differences between the two models. In 
Panel (a), both the reference and test models have the same ra- 
dius, 695.99Mm, the difference between the two models is the 
equation of state used. In Panel (b), the models have identical 
input physics, but radii of the two models are different. The ref- 
erence model has a radius of 695.99 Mm while the test model has 
a radius of 695.78 Mm. 

f-mode frequencies. In Sections 3 and 4 we discuss the effect 
of the radius change on the deduced depth of the solar con- 
vection zone and helium abundance respectively. In Section 
5 we try to investigate if an error in equation of state can 
explain the remaining difference between the sound speed 
of the Sun and solar models inside the convection zone. We 
present the conclusions from this study in Section 6. 



2 INVERSION RESULTS 

In this work we have used the subtractive optimally lo- 
calized averages (SOLA) method of Pijpers & Thompson 
(1992) to invert Eq. (1). The procedure as applied to de- 
termine the sound-speed difference is discussed in Basu et 
al. (1996). 

Three sets of helioseismic data have been used for this 
work. The first set is the data from the Big Bear Solar Obser- 
vatory (Libbrecht, Woodard and Kaufman, 1990) combined 
with the low degree (0 < I < 3) modes from the Birmingham 
Solar Oscillation Network (Elsworth et al. 1994), this set is 
referred to as the BBSO set. The second set consists of the 
frequencies obtained from the average spectra of months 4 - 
14 obtained by the Global Oscillation Network Group (here- 
after referred to as the GONG data), and the third set are 
the frequencies from the data obtained by the SOI/MDI in- 
strument on board the SOHO spacecraft during its first 144 
days in operation (Rhodes et al. 1997), which we shall refer 
to as the MDI data. We restrict ourselves to the frequency 
range of 1 to 3.5 mHz. The BBSO set is further restricted to 
a degree range of I < 140 to avoid frequencies determined 
by the ridge fitting method which may have large systematic 
errors. 

We have used two reference models with identical physi- 
cal inputs, but different radii for the inversions. Both models 
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Fig. 2 The relative sound-speed difference between the Sun and 
two reference models R99 and R78 at fixed fractional radius. The 
differences are in the sense (Sun— modcl)/model. The triangles are 
the results obtained with a reference model of R99 (radius 695.99 
Mm), the circles are obtained with reference model R78 (radius 
695.78 Mm). The vertical error bars show the la error in the 
inversions due to errors in the data. The horizontal error-bars are 
a measure of the resolution. The sound-speed difference between 
the two models can be seen in Fig. 1(b). 

are static-models, constructed with the OPAL equation of 
state (Rogers, Swenson & Iglesias 1996) and OPAL opacities 
(Iglesias & Rogers 1996) which are supplemented by the low 
temperature opacities of Kurucz (1991). The CM formula- 
tion is used to calculate the convective flux. To construct 
the model we use the helioseismically determined hydrogen 
abundance profile from Antia & Chitre (1997a), and also 
use the heavy elements abundance profile adopted by them. 
The surface ratio of heavy elements to hydrogen abundance 
was constrained to be 0.0245, the value determined from 
observations by Grevesse & Noels (1993). The radii used 
are the standard value of 695.99 Mm (Model R99) and the 
reduced value 695.78 Mm (Model R78) as determined by 
Antia (1997). In this paper we denote the radius of the solar 
models as i? s , to avoid confusion with Rq, which is normally 
assumed to have the fixed value of 695.99 Mm. 

The sound-speed difference between the Sun and the 
two models are shown in Fig. 2. The sound-speed difference 
between the two models has been shown earlier in Fig. 1(b), 
and hence it is easy to see that there is a real difference be- 
tween solar sound-speed deduced by using the two models. 
The same is true for solar density. To make comparisons eas- 
ier, we have represented the calculated sound-speed profiles 
as differences with respect to the Model S of Christensen- 
Dalsgaard et al. (1996) in Fig. 3(a). 

Thus it can be seen that there is a significant difference 
between the results obtained using the two models, which 
arises because of the difference in radius. In fact, for sound- 
speed, the relative difference between the two solutions is 
of the same order as the relative sound-speed difference be- 
tween the Sun and the two reference models. Thus it is im- 
portant to invert solar oscillation frequencies using reference 
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Fig. 3 (a) The sound-speed and (b) density inside the Sun as 
obtained by the two reference models R99 and R78 shown with 
respect to Model S of Christensen-Dalsgaard et al. (1996). In 
both panels, the continuous line represents the result obtained 
with model R78 (R s = 695.78 Mm) and the dashed line that that 
obtained using model R99 (R B = 695.99 Mm). The dotted lines 
represent the la error envelope. 

models with the same radius as the Sun. Fig. 3(b) shows the 
density results obtained using the two models. Again there 
is a real difference in the result obtained, though in the case 
of density, the changes are within the la error limits. 

The sound-speed difference near the surface between 
Model S and the Sun shows a very sharp dip which may 
be attributed to the difference in radius. In fact, the sound- 
speed difference at fixed radius between models R99 and 
R78 shows a sharp decrease at the surface, which may lead 
one to assume that the change in the radius of the refer- 
ence model will reduce the dip. However, this decrease is 
not very evident in differences taken at fixed fractional ra- 
dius (cf., Fig. lb), which is what is used in the inversions, 
and hence it is unlikely to affect the inversion results signifi- 
cantly. In fact, we can see from Figs. 2 and 3 that the surface 
behaviour of the solutions is not very different regardless of 
the radius of the model used for inverting the solar oscilla- 
tion frequencies. Thus we believe that the sharp decline in 
the relative sound-speed difference is more a result of the 
other physical properties of the model. 

To confirm this point we have constructed reference so- 
lar models using physical inputs which give rise to large 
changes in the surface layers as compared to the two refer- 
ence models used above. One set consists of models with 
the two different radii and constructed using the mixing 
length theory (MLT) for calculating the convective flux in- 
stead of CM formulation. The models are called MLT99 (for 
R s = 695.99 Mm) and MLT78 (for R s = 695.78 Mm). We 
also have one low radius model where instead of Kurucz 
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Fig. 4 (a) The relative sound-speed difference between the refer- 
ence models R99 and R78 and the models with modified physics. 
The differences are calculated in the sense (reference — modified) . 
The continuous line is the difference between models R78 and 
OP78, the dotted line is the difference between models R99 and 
MLT99 and the dashed line is the difference between models R78 
and MLT78. Note that the dotted and dashed lines more or less 
coincide, (b) The relative sound-speed difference between the Sun 
and model R99 (triangles) and MLT99 (squares), (c) The relative 
sound-speed difference between the Sun and models R78 (trian- 
gles), MLT78 (squares) and OP78 (circles). The error-bars are 
omitted in this panel for the sake of clarity. 

opacities, the OPAL opacities as used as far as the table 
continue. Below that opacities from Alexander (1975) are 
used. However, in that region the temperature stratification 
is obtained using the T — t relation of Vernazza et al. (1981). 
This model is called OP78. 

The relative difference in sound-speed between our orig- 
inal reference models and the modified models are shown 
in Fig. 4(a). Note that the modified models have a higher 
sound speed in the surface layers as compared to the models 
with the CM formulation of convection and Kurucz opaci- 
ties, thus the sound-speed difference between the Sun and 
the modified models is likely to show a dip at the surface. 
The results of the inversion are shown in Fig. 4 (b) and (c). 
The new models, irrespective of radius, do show a sharp de- 
crease in the sound-speed difference in the surface layers. 
Thus it is clear that a difference in radius is not the cause 
of the near-surface plunge in the relative sound-speed differ- 
ence seen between the Sun and Model S. 

Fig. 4 also supports the result of Basu & Antia (1994a) 
that the models constructed with the CM formulation are 
closer to the Sun than those with MLT formulation. It also 
appears that the Kurucz (1991) opacities fit solar data better 
than just the OPAL opacities supplemented by those from 
Alexander (1975). This appears to support the results of 



Antia & Basu (1997) using just the frequency differences. 
However, there remains a small dip near the surface which 
points to remaining deficiencies in model R78. 

The change in the deduced solar sound-speed pro- 
files because of changes in the radius necessitates the re- 
calibration of quantities derived from the sound-speed dif- 
ference. We therefore, look next at the depth and helium 
abundance of the solar convection zone. 



3 DEPTH OF THE CONVECTION ZONE 

The transition of the temperature gradient from the adia- 
batic to radiative values at the base of the solar convection 
zone (CZ) leaves its signature on the sound speed. Thus 
helioseismic measurements of the sound speed enables a de- 
termination of the position of the base of the convection 
zone (cf., Christensen-Dalsgaard, Cough & Thompson 1991; 
Kosovichev & Fedorova 1991; Guzik & Cox 1993; Basu & 
Antia 1997). 

3.1 The Procedure 

We use the technique described in Basu & Antia (1997) to 
determine the position of the CZ base, 7V To recapitulate 
briefly, if there are two otherwise similar solar models with 
different depths of the convection zone, then just below the 
base of the convection zone the model with a deeper con- 
vection zone will have a larger sound speed than the other. 
This observable difference of sound speeds can be calibrated 
to find the convection-zone depth of a test model or the Sun. 
Asymptotically, the frequency differences between a so- 
lar model and the Sun, or between two solar models can 
be written as (Christensen-Dalsgaard, Cough & Thompson 
1989) 

S(w)— =H 1 {w) +H 2 {u>), (2) 

where w — u>/(£ + 0.5), and S(w) is a known function for 
a given reference model. The functions Hi(w) and H^iui) 
can be found by a least-squares fit to the known frequency 
differences. Hi(w) can be inverted to obtain the sound-speed 
difference, <5c/c, between the reference model and the Sun. 
However, that is not required as Hi (w) itself can be used to 
determine the convection-zone depth. 

If (j>( w ) i s the -Hi(w) between two solar models which 
differ only in the depth of the convection zone, then H\ (w) 
for any other pair of models can be written as 



Hi(w) = P<f>(w) +H s (w) 



(3) 



where H 3 (w) is a smooth component of Hi (w) which results 
from sound-speed differences that arise from differences in 
the equation of state, abundances, surface physics etc., and 
the first term is the contribution to Hi (w) due to the sound- 
speed difference caused by the difference in rt, the position 
of the base of the CZ. Thus if j3 is determined by a least 
squares fit, the unknown rb of the Sun can be obtained. 

We determine f3 for a series of calibration models with 
different n,, and interpolate to find the position for which 
/3 = 0. The error in the CZ depth arising from those in the 
observed frequencies is determined by Monte Carlo simula- 
tions. 
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As is clear from the above discussion, we need to use 
models with a specified depth of CZ. This is easier with mod- 
els of the solar envelope, where position of the CZ base and 
the helium abundance of the solar envelope can be specified. 
As Basu & Antia (1997) have shown, the largest source of 
systematic error in the results is the hydrogen abundance 
profile. This is because the sound speed near the base of 
the CZ is affected not only by the change in the tempera- 
ture gradient, but also by the change in the mean molecular 
weight due to gravitational settling of helium. Since the X 
and Z profiles are not known exactly, we use models with 
three types of profiles: 
ND: The hydrogen and heavy element abundance profiles 
from the no-diffusion model of Bahcall & Pinsonncault 
(1992). 
DIF: The hydrogen abundance profile from solar model of 
Bahcall & Pinsonncault (1992) incorporating the dif- 
fusion of helium and Z profile from Proffitt (1994). 
INV: The hydrogen abundance profile obtained from helio- 
seismic inversions (Antia & Chitre 1997a) and the Z 
profile adopted in that paper. 
For each type of model, we have two sets of calibra- 
tion models, one with the standard radius of 695.99 Mm 
and the other with the reduced radius of 695.78 Mm. Each 
calibration set consists of 5 models with CZ-base positions 
at 0.709, 0.711, 0.713, 0.715 and 0.717 R s . Thus we have 6 
sets of calibration models. In addition, we have test mod- 
els for each hydrogen abundance profile and each radius, 
with Tb = 0.712 R s . We also have test models with an even 
lower radius of 695.57 Mm, to check for trends. For conve- 
nience, the models are referred to by their radius (99 for 
695.99 Mm, 78 for 695.78 Mm and 57 for 695.57 Mm), hy- 
drogen and heavy element abundance profile and radius of 
the CZ base, e.g., DIF99-0.709 is the model with the DIF 
X and Z profiles, with rt = 0.709i? s and standard radius 
(R s = 695.99 Mm), INV78-0.711 s the model with INV X 
profile, r;, = 0.711i?s and smaller radius (Rs — 695.78 Mm), 
etc. All models have been constructed with the OPAL equa- 
tion of state, OPAL and Kurucz opacities and with the CM 
formulation for calculating convective flux. They have iden- 
tical surface X of 0.7342 and surface (Z/X) of 0.0245. The 
systematic errors arising due to variations in these parame- 
ters were studied by Basu and Antia (1997). 

3.2 Results 

The results for the test models are listed in Table 1. We 
have shown the results for the INV models only. The results 
for the other calibration models are similar. We can see that 
there is indeed some difference in the measured value of ra- 
dial position of the CZ base, when reference models with 
a different radius are used. The difference is of the order of 
the relative difference in the radius, but is somewhat smaller 
than error in the sound speed caused by mismatch of radius. 
The reason for this can be seen in Fig. 5, where we have 
shown the function H\(w) for model INV78-0.712 and the 
INV99 and INV78 calibration models. Hi(w) is flat in the 
CZ for INV78 calibration models. This is consistent with 
the convection-zone sound-speed difference between these 
models and the test model INV78-0.712 being 0. The pre- 
dominant difference between the results with INV99 and 
INV78 models is a linear shift. The truly linear part of the 



Table 1. CZ results for the test models 



Test model 


Deduced r^/Rs 




Rs = 695.99Mm Rs = 695.78Mm 


Calib. models with INV X profile 


INV99-0.712 


0.712000 0.712258 


INV78-0.712 


0.711746 0.712000 


INV57-0.712 


0.711498 0.711746 
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Fig. 5 The function Hi(w) between the model INV78-0.712 and 
(a) the INV99 calibration models and (b) the INV78 calibration 
models. In both panels, the continuous line shows Hi(w) between 
the test model and the calibration model with r\, = 0.709i? s , 
the dotted line is with r^ = 0.711-R S , small-dashed line with 
r;, = 0.713-Rs, long-dashed line with rj, = 0.715-R S and dot-dashed 
line with r b = 0.717R B - 




Fig. 6 The function Hl(w) between the Sun and (a) the INV99 
calibration models and (b) the INV78 calibration models. The 
line types are the same as in Fig. 5. 

shift is easily removed by the function H s (w) in the fit to 
Eq. (3), hence the r b results for this model are not very dif- 
ferent. Nevertheless, the error arising due to difference in 
radius is larger than the error in the solar results because 
of errors in the frequencies, and hence the solar CZ depth 
should be determined using the models with the correct ra- 
dius. The difference in rt, due to a difference in the radius is 
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Table 2. CZ results for the Sun 



Data set 



Deduced r b /R s 



ND78 models 



INV78 models 



DIF78 models 



BBSO 
GONG 

MDI 

Weighted Average 



0.714219 ±0.000179 
0.714149 ± 0.000062 
0.714443 ± 0.000080 
0.714267 ± 0.000050 



0.713435 ± 0.000175 
0.713430 ± 0.000060 
0.713601 ± 0.000078 
0.713494 ± 0.000049 



0.710665 ±0.000181 
0.710631 ± 0.000062 
0.710868 ± 0.000081 
0.710722 ± 0.000051 



much smaller than that due to differences in the hydrogen 
abundance profile. For example, the error in determining the 
position of the CZ base can be as much as 0.0035 Rs (cf., 
Basu & Antia 1997) between the DIF and ND profiles. 

The function Hi(w) between the Sun and the two sets 
of INV calibration models are shown in Fig. 6. Note that 
Hi(w) for models with R s — 695.78 Mm is much flatter in 
the CZ as compared to models with R s = 695.99 Mm. This 
gives added confidence that the solar radius is indeed close 
to 695.78 Mm. Table 2 lists the position of the solar CZ base 
as obtained using different data sets and different calibration 
models. 

Since it can be shown that models with the DIF pro- 
file have too sharp a gradient in the X-profile at the CZ 
base (Basu 1997a), we can disregard the value obtained with 
those models. There is an enormous amount of evidence that 
diffusion of helium and heavy elements does take place below 
the solar CZ base (see e.g. Basu 1997b), the value obtained 
with the ND models is thus an upper limit to the value of 
r b , as models with diffusion will give lower values. The value 
obtained with the INV models is probably the most reli- 
able determination. The change in the results is negligible if 
instead of the INV profiles, those obtained from stellar evo- 
lution codes that incorporate mixing due to rotation (e.g., 
Richard et al. 1996) are used. 

The various other systematic errors that can creep into 
the determination of the position of the CZ base have been 
discussed in detail in Basu & Antia (1997). They are in- 
variably larger than the statistical error due to the error in 
frequencies. This gives an error of 0.0005_R S , which is larger 
than the error due to the uncertainty in the radius. 

The base of the convection zone is therefore estimated to 
be at 0.7135 ± 0.00057? s . If models with the standard radius 
are used, we find a position of 0.7132 ± 0.00057? s using the 
same data sets. 

3.3 CZ depth and Z/X from density differences 

The sound speed inside the convection zone is not partic- 
ularly sensitive to the composition or the position of the 
convection-zone base. The difference in the sound speed be- 
tween various models occurs below the CZ base. Thus it is 
difficult to disentangle the change in sound-speed below the 
CZ base because of change in temperature gradient from 
that due to changes in the mean-molecular weight due to 
gravitational settling of helium and heavy elements. This is 
the reason that the models with different composition pro- 
files give seemingly different estimates of the position of the 
CZ base. The density in the CZ is, however, affected by the 
position of the CZ base, which arises because of difference 
in mixing length required to achieve different depths of CZ. 
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Fig. 7 The relative density difference between the Sun and the 
INV78 calibration models. The thin-continuous, small-dashed, 
long-dashed, dot-small dashed and dot-long dashed lines are the 
differences for the models with r b = 0.709, 0.711, 0.713, 0.715 and 
0.717 R s respectively. The thick continuous line is for the model 
with the CZ base at the position estimated using the above cal- 
ibration models (r b = 0.7135R S ). The dotted lines show the la- 
error in the density inversion. 



This can be seen from Fig. 7, where we show the density 
difference between Sun and the INV calibration models. We 
have also shown the density difference for the model with 
the CZ base at the deduced position of 0.7135i? s . Note that 
the density difference for this model is consistent with 0, 
which gives us added confidence in the results. It may be 
noted that the position of the base of the convection zone is 
determined by the opacity and in principle, it is possible to 
estimate the opacity at the CZ base by using the inversion 
results coupled with an estimate of the position of the CZ 
base. Analyses using density differences unfortunately suffer 
from the drawback that density in the CZ is also dependent 
on the low-temperature opacities used and to some extent 
on the composition within the CZ. 

As can be seen from Fig. 7, the density difference be- 
tween the INV78 model with r b — 0.7135 and surface 
Z/X = 0.0245 and the Sun is almost zero. This would mean 
that assuming that r b and Z/X are correct, the opacities 
used in the model, i.e., OPAL supplemented by Kurucz, arc 
consistent with helioseismic data. The uncertainty in the es- 
timated depth of the CZ base i.e., 0.0005i? s corresponds to 
a change of opacity at the CZ base by 0.6%, which is much 
smaller than the probable uncertainties in the opacity. The 
largest change in opacity occurs when Z/X is changed. If 
Z/X is changed by the quoted uncertainty of 10%, the opac- 
ities also show almost 10% change. Thus we can conclude 
that with the current values of Z/X the OPAL opacity ta- 
bles are consistent with the estimated opacity at the base of 
the solar CZ. 
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The good agreement between the estimated opacity and 
that calculated from the OPAL tables using the currently ac- 
cepted value of Z/X suggests that the uncertainty in Z/X 
should be much smaller than the quoted error of 10%. As- 
suming the correctness of the OPAL opacities we can use 
the CZ-density differences between the Sun and the models 
to put limits on the amount of heavy elements in the solar 
envelope. The solar models constructed with rt, = 0.7135i? s , 
Y = 0.2478 and Z/X = 0.0245 shows a near perfect fit to 
the CZ density. If we keep Y the same, then we find that a 
change in Z/X of ±0.0001 changes the density difference by 
±1(7. If we keep X fixed and change Z/X then a change of 
±0.0002 changes density by ±lo\ In addition there are un- 
certainties caused by errors in the depth of the CZ, the low 
temperature opacities, the hydrogen abundance in the CZ, 
and uncertainties in the solar-radius. Taking into account 
all these uncertainties we estimate that the heavy element 
abundance in the solar envelope is Z/X — 0.0245 ± 0.0006, 
i.e., the uncertainty is much less than the normally assumed 
uncertainty of 10% in spectroscopically estimated value. 



3.4 Overshoot below the solar CZ base 

Solar oscillation frequencies can be used to determine 
the extent of overshoot below the solar convection zone if it 
is assumed that the overshoot region is adiabatically strati- 
fied and that there is an abrupt transition to radiative strat- 
ification below the overshoot zone. This change in stratifi- 
cation introduces a discontinuity in the first derivative of 
the sound-speed and leaves a characteristic oscillatory sig- 
nal in the frequencies. The amplitude of the signal can be 
calibrated to determine the extent of overshoot (cf., Mon- 
teiro, Christensen-Dalsgaard & Thompson 1994; Basu & An- 
tia 1994b; Basu 1997a). Basu (1997a) had found an upper 
limit of 0.05H P on the overshoot region. That was done us- 
ing models with the standard radius of 695.99 Mm. We have 
tried to see whether models with the lower radius show a 
change in the amplitude of the signal and whether the use 
of the new data changes the limits. The procedure used for 
fitting is described in Basu (1997a). 

We constructed a few models to test the change of the 
amplitude with radius and extent of overshoot. All mod- 
els are INV models. The CZ depth of the models is ap- 
propriate for the radius of the model. The results for the 
different models and data sets are shown in Table 3. The 
addition of the two data sets does not change the observed 
results significantly. The change in the amplitude of the sig- 
nal with change in radius is smaller than the statistical error 
in the amplitude due to uncertainties in observed frequen- 
cies. Hence the conclusions about the extent of overshoot 
below the solar convection zone are not affected by possible 
error in solar radius. 



4 HELIUM ABUNDANCE 

Spectroscopic measurements of the abundance of helium in 
the Sun being very uncertain, helioseismology plays a major 
role in determining the helium abundance in the solar enve- 
lope. The abundance is obtained from the variation of the 
adiabatic index of the solar material in the helium ionisation 
zone, which in turn affects the sound speed. 



Table 3. Amplitude of signal form CZ base 



Model 


Amplitude (^Hz) 
INV99 INV78 


Ov.=0 H p 
Ov.=0.05 Hp 
Ov.=0.10 H p 


0.8020 
0.9014 

1.0428 


0.8073 
0.9071 
1.0494 


Observations 



BBSO 

GONG4-14 

MDI 

Weighted Average 



0.879 ±0.106 
0.782 ±0.032 
0.771 ±0.034 
0.790 ± 0.025 



4.1 The procedure 

The technique used here is the same as that in Antia & 
Basu (1994b) and Basu & Antia (1995), and is very similar 
to that used to determine the depth of the CZ zone. The 
starting point of the procedure is again Eq. (2). However, 
since the data do not extend to high degrees and the data 
sets have hardly any modes with turning points in the helium 
ionisation zone (r > 0.98i? s ), it is difficult to use Hi(w) 
to determine the helium abundance. We have to use the 
function H^iuS) instead. It has been shown that (e.g., Basu 
& Antia 1995; Perez Hernandez & Christensen-Dalsgaard 
1994), the helium ionisation zone being sufficiently close to 
the solar surface leaves its imprint on the function H^iui) 
also, in the form of a hump around 2.5 mHz, which can be 
calibrated to determine the helium abundance. The effect of 
difference in other physical properties, like the equation of 
state, manifests itself as a smooth curve on which the hump 
due to difference in helium abundance is superimposed. 

Thus, analogous to Eq. (3), we write the function H^iuj) 
between the Sun and a calibration model of known helium 
abundance as: 



ff 2 (w)=/?$(w) + # s (w), 



(4) 



where $(w) is the difference in H^iuf) between two models 
which differ only in the abundance of helium. 

For this work too it is preferable to use envelope models 
since the calibration models are required to have specified 
values of the helium abundance. Since by far the largest 
error in helium abundance determination occurs due to a 
mismatch of the equation of state (cf., Antia & Basu 1994b), 
we use models with both OPAL and the MHD (Hummer & 
Mihalas 1988; Mihalas, Dappen & Hummer 1988; Dappen 
et al. 1988) equations of state. For each equation of state, 
we have constructed calibration sets with radii 695.99 Mm 
and 695.78 Mm. Each calibration set consists of five models 
with hydrogen abundance X of 0.68, 0.70, 0.72, 0.74 and 
0.76. The value of the metal abundance Z was kept fixed 
at 0.018 for the OPAL models. For the MHD models, the 
available tables force Z — 0.02. All models have the CZ base 
at Tb — 0.713 R s . In addition, we have test models with 
X — 0.73 for both equations of state and both values of 
the radius. We also have constructed test models with R s = 
695.57 Mm. The models are labelled by their EOS (OPAL 
or MHD), radius (99 for 695.99 Mm, 78 for 695.78 Mm and 
57 for 695.57 Mm), and surface hydrogen abundance. Since 
the helium ionisation zone is very far from the CZ base, it is 
not expected to be affected by the abundance profiles below 
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Table 4. Helium abundance results for the test models 



Test model 



Exactr 



a 



Deduced Y 
695.99 Mm _R S = 695.78 Mm 



Calib. models with OPAL EOS 



OPAL99-0.73 
OPAL78-0.73 
OPAL57-0.73 
MHD78-0.73 



0.252 
0.252 
0.252 
0.250 



0.25184 
0.25212 
0.25240 
0.24938 



0.25156 
0.25184 
0.25212 
0.24904 



Calib. models with MHD EOS 



MHD99-0.73 
MHD78-0.73 
MHD57-0.73 
OPAL78-0.73 



0.250 
0.250 
0.250 
0.252 



0.24987 
0.25016 
0.25046 
0.25373 



0.24957 
0.24987 
0.25016 
0.25355 



Table 5. Helium abundance results for the Sun 


Data set 


Deduced Y 
R s = 695.99 Mm R s = 695.78 Mm 


OPAL models 


BBSO 
GONG 

MDI 


0.2476 ± 0.0002 0.2474 ± 0.0002 
0.2471 ± 0.0001 0.2470 ± 0.0001 
0.2488 ± 0.0001 0.2488 ± 0.0001 


MHD models 


BBSO 
GONG 

MDI 


0.2514 ±0.0002 0.2510 ±0.0002 
0.2499 ± 0.0001 0.2496 ± 0.0001 
0.2524 ± 0.0001 0.2522 ± 0.0001 



the CZ base. Hence, the models constructed for calibrating 
the helium abundance do not assume any diffusion of helium 
or heavy elements below the CZ base. 

4.2 Results 

The value of Y obtained for different models are listed in Ta- 
ble 4. Note that unlike the case of the CZ depth, where the 
CZ depth of the models with same i? s could be determined 
accurately, there is an error of about 0.0002 in the determi- 
nation of the helium abundance even when the radius is the 
same. The main reason is that the signature due to differ- 
ences in helium abundance is not as clear-cut as that due to 
differences in the position of the CZ base. Another reason 
for this is that the mode-sets used do not have modes with 
turning points in the helium ionisation zone. Nonetheless, 
we see that there is an additional error if the radius of the 
calibration and test models differs. The magnitude of the 
difference is of the order of the fitting errors in the deter- 
mination of helium abundance. However, as we shall show 
below, there are larger errors caused by a number of other 
factors, such as opacity used etc. Hence the solar results are 
unlikely to be affected significantly by a change in radius of 
the calibration models. 

Since the contribution to H2(uj) at any given frequency 
comes from modes with different turning points, the deduced 
helium abundance may be affected by model differences not 
necessarily localised at the helium ionisation zone. To test 
this we have constructed a few additional models to test the 
error in the deduced helium abundance caused because of 
differences in composition profile below the CZ base, depth 
of CZ and opacities used. As expected, the composition pro- 
file below the CZ base does not affect the determination of 
He abundance. However, the depth of the CZ does. An error 



of 0.0005i? s in the CZ depth causes an error of 0.0003 in Y. 
A much larger error of 0.0006 is produced if the low tem- 
perature Kurucz (1991) opacities are replaced by those of 
Alexander (1975). We estimate that the uncertainty in the 
deduced helium abundance due to all these factors is about 
0.001. 

The results for the Sun are listed in Table 5. The error- 
estimate corresponds merely to the error in the result due 
to those in the frequencies, and as can be seen in light of the 
discussion above, are much smaller than the systematic er- 
rors involved. Unlike in the case of CZ depth determination, 
the different data sets give results that are not consistent 
with each other within the errors. We believe this is a re- 
sult of the different degree ranges of the three data sets. 
For a definitive determination of the helium abundance, one 
needs high-degree p-modes which have lower turning points 
in the helium ionisation zones. None of the data-sets used 
has such high degree modes. The highest-degree p-mode in 
the MDI set has I = 194, in the GONG set it has £ = 150 
and the BBSO set has been deliberately restricted to modes 
of I < 140 to avoid using modes with large systematic er- 
rors. Since none of the mode-sets has modes which sample 
the helium ionisation zone very well, it is not very surprising 
that the results with the three data sets are different 

As expected (cf., Basu & Antia 1995; Kosovichev 1996; 
Baturin & Ayukov 1997), MHD and OPAL models give dif- 
ferent results for the Sun. The result obtained for the OPAL 
models, an average of the three sets being Y = 0.2478, is not 
very different from the earlier result of Y — 0.249 obtained 
by Basu & Antia (1995), and are the same as those of Koso- 
vichev (1996). The current results are, however, expected 
to be more reliable since the radius of the models is consis- 
tent with that of the Sun and the data-sets used are better. 
For models with the standard radius of 695.99 Mm, we find 
Y — 0.2479, with an uncertainty of 0.001 in the results due 
to fitting errors and systematic error mentioned above. 

The result for the MHD models on the other hand is 
somewhat higher than that obtained by Basu & Antia (1995) 
but are consistent with those of Baturin & Ayukov (1997). 
The results obtained by different workers using MHD mod- 
els have a wide dispersion. This may be due to the fact that 
the sound-speed difference between all MHD models and 
the Sun show a distinct deviation just below the helium ion- 
isation zone (cf., Dziembowski, Pamyatnykh & Sienkiewicz 
1992; Basu & Antia 1995), which can contaminate the signal 
from the helium abundance to varying degrees depending on 
the procedure followed to estimate the helium abundance in 
the solar envelope. 



5 
CZ 



THE SOUND-SPEED DIFFERENCE IN THE 



Fig. 8 shows the sound-speed difference between the 
Sun and four solar models in the convection zone. There 
are two OPAL models and two with the MHD EOS. For 
each EOS, one model has the standard radius of 695.99 Mm 
and one the reduced radius of 695.78 Mm. Each model has 
a CZ-base position and helium abundance appropriate to 
it's radius and EOS (cf., Sections 3 and 4). Note that the 
sound-speed of the Sun is higher than that of the two models 
with the standard radius. For the reduced-radius models, 
the solar sound-speed is higher than the sound-speed of the 
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Fig. 8 The relative sound-speed difference between solar envelope 
models and the Sun. The continuous line and the small-dashed 
line are for models with R a = 695.78 Mm, constructed with OPAL 
and MHD equations of state respectively. The dotted lines are the 
ler error envelope. The long-dashed and the dot-dashed lines are 
for models with R B = 695.99 Mm, constructed with OPAL and 
MHD equations of state respectively. For the sake of clarity, the 
error envelopes have not been shown for these models. 

MHD model, but is marginally lower than that of the OPAL 
model. The overall difference is however much lower than the 
models with the standard radius. 

There are a number of possibilities as to why the sound- 
speed in the lower part of the CZ in the models do not match 
that of the Sun, these are, for instance, that the change in 
radius has not been estimated correctly or that the equation 
of state is not correct. 

Looking at Fig. 8, it appears that the change in radius 
has over-corrected the sound-speed in the case of the OPAL 
model, while the correction is not enough for the MHD mod- 
els. Thus, for the OPAL models the reduction of radius to 
695.78 Mm is an overestimate of the reduction, while for 
the MHD models it is an underestimate. The radii at which 
the models show very little difference with respect to the so- 
lar sound speed in the lower CZ are 695.82 Mm and 695.67 
Mm respectively for the OPAL and MHD EOS. This addi- 
tional change in radius, however, causes the f-mode frequen- 
cies of the model to deviate from the observed frequencies. 
This means that either the surface layers of the solar mod- 
els have not been modelled properly, or other processes, like 
flows etc., may need to be postulated to explain the f-mode 
frequency difference. 

The structure of lower part of the solar convection zone 
is determined mainly by the equation of state, hence an er- 
ror in the equation of state will give rise to a sound-speed 
difference between a solar model and the Sun. Thus if we 
assume that the radius and helium abundance in CZ have 
been estimated correctly, the remaining discrepancy in the 
sound speed may be due to possible error in the equation of 
state. It is clear from Fig. 9 that the OPAL EOS is much 
closer to that of the solar material than the MHD EOS. How- 
ever, there is still some remaining difference. We find that 
for the reduced radius models, a large part of the difference 
can be eliminated if the partial derivative (dp/dT) p for the 
OPAL equation of state is lowered by a factor (1 — 0.00027), 
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Fig. 9Thc relative sound-speed difference between solar envelope 
models and the Sun. Both models have a radius of 695.78 Mm. 
The continuous line is for a model constructed with a modified 
OPAL EOS, and the dashed line is for a model constructed with 
modified MHD EOS. 

and the other thermodynamic quantities like C p and Li are 
changed consistently. This corresponds to a relative differ- 
ence of about —0.00022 in Ti between the two models in 
the lower part of the CZ. For the MHD EOS, the partial- 
derivative should be larger by a factor 1.00053, this results 
in a relative difference of 0.00043 in Ti. This small differ- 
ence largely eliminates the sound-speed difference in the CZ, 
as can be seen from Fig. 9. It is not obvious that the dis- 
crepancy is indeed due to an error in equation of state, but 
these differences probably give an estimate of uncertainties 
in current equations of state. This change in EOS does not 
change the f-mode frequencies significantly. 



6 CONCLUSIONS 

We have shown that a mismatch in the radius of the Sun and 
that of reference solar models can introduce considerable 
error in the inverted sound-speed and density profiles of the 
Sun. The relative difference in the sound speed obtained by 
inversion due to the change in radius is of the same order 
as the relative difference between the sound speed of the 
Sun and the solar models used. Therefore, inversion of solar 
oscillation frequencies to determine solar structure should be 
performed using models with radius as close as possible to 
the correct value. Thus it is necessary to obtain an accurate 
estimate of the solar radius. 

We find that regardless of the radius of the solar models 
used, solar models constructed using the CM formulation of 
convection give a better fit to the near-surface solar sound 
speed than models constructed with the mixing length for- 
malism. Models constructed with the mixing length formal- 
ism show a sharp dip in the relative sound-speed difference 
with respect to the Sun in the outer regions, and this dip is 
not removed by the reducing the radius of the solar models. 

We have used new data-sets to check how the change 
in radius results in a small change in the measured depth of 
the solar convection zone and in the helium abundance. We 
obtain a radial distance of (0.7135 ± 0.0005)i? s for the base 
of the CZ. We also find that the OPAL tables are consistent 
with the solar opacities at the CZ base for a surface Z/X of 
0.0245. Assuming that the OPAL opacities are correct, wc 
deduce that the the surface Z/X = 0.0245 ± 0.0006 with an 
error estimate which is much less than the usually quoted 
value of 10%. The limit on the extent of overshoot does not 
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change because of the small change in radius, and remains 
at 0.05H P (2800 Km). 

Using the new data sets and models constructed with 
the OPAL equation of state, we find that the helium abun- 
dance in the solar envelope is 0.2478±0.001. Once again, the 
estimate of helium abundance is not significantly affected by 
the error in radius as the effect is smaller than the systematic 
errors in this measurement. 

If the correct R s is indeed 695.78 Mm and abundances 
in CZ have been estimated correctly, then the remaining dis- 
crepancies between the CZ sound-speed profile of the models 
and the Sun can be explained as errors in the equation of 
state used. It can be shown that by changing (dp/dT) p for 
each EOS suitably, one can reduce the difference. For OPAL 
models this has to be reduced by a factor of 0.99973, while 
for MHD models this has to be increased by a factor 1.00053. 
Changes in (dp/dT) p automatically changes the adiabatic 
index Pi, and corresponds to a decrease in Pj by a factor 
of 0.99978 for OPAL models and an increase by a factor of 
1.00043 for MHD models. This probably gives an estimate 
of uncertainties in the current equations of state. 
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